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Constituents of complex systems interact with each other and self-organize to 
form complex networks. Empirical results show that the link formation pro¬ 
cess of many real networks follows either the global principle of popularity 
or the local principle of similarity or a tradeoff between the two. In partic¬ 
ular, it has been shown that in social networks individuals exhibit significant 
homophily when choosing their collaborators. We demonstrate, however, that 
in populations in which there is a division of labor, skill complementarity is 
an important factor in the formation of socioeconomic networks and an in¬ 
dividual’s choice of collaborators is strongly affected by heterophily. We an¬ 
alyze 124 evolving virtual worlds of a popular “massively multiplayer online 
role-playing game” (MMORPG) in which people belong to three different pro¬ 
fessions and are allowed to work and interact with each other in a somewhat 
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realistic manner. We find evidence of heterophily in the formation of collabo¬ 
ration networks, where people prefer to forge social ties with people who have 
professions different from their own. We then construct an economic model 
to quantify the heterophily by assuming that individuals in socioeconomic sys¬ 
tems choose collaborators that are of maximum utility. The results of model 
calibration confirm the presence of heterophily. Both empirical analysis and 
model calibration show that the heterophilous feature is persistent along the 
evolution of virtual worlds. We also find that the degree of complementarity in 
virtual societies is positively correlated with their economic output. Our work 
sheds new light on the scientific research utility of virtual worlds for studying 
human behaviors in complex socioeconomic systems. 


Introduction 

Complexity emerges in the evolving and self-organizing processes of many natural, social, tech¬ 
nological, and biological systems. The constituents of a complex system interact with each other 
and form complex evolving networks, where the constituents are nodes and their interaction re¬ 
lationships are links ( 1 - 6 ). For many real networks, the link formation process follows either 
the global principle of popularity in which a node tends to link with high-degree nodes (7,5), or 
the local principle of similarity in which a node tends to link with nodes having traits similar to 
its own (9), or a tradeoff between them (9). 

In the sociological literature the local principle of similarity, i.e., the phenomenon that “birds 
of a feather flock together,” is known as homophily ( 10 ). There is much empirical evidence in¬ 
dicating that individuals prefer to forge social ties with people whose traits such as education, 
race, age, and sex are the same as their own ( 11 - 14 ). Such homophilous behaviors are ubiq¬ 
uitous in social networks and have been well documented ( 10 - 12 , 14 - 18 ). In addition, the 
similarity shared by individuals in a group is often a significant predictor of a group’s altruism 
level and its ability to cooperate ( 19 ). Sociological literature argues that human societies tend to 
display two social systems: (i) homophilous, in which people seek out people who are similar, 
and (ii) heterophilous, in which people seek out people who are different ( 20 ). The evidence 
indicating the actual existence of heterophilous societies is rare, however. 

In general, it has long been accepted that one of the most significant factors in increasing 
productivity in modem human societies has been the division of labor ( 21 ). Thus we might 
assume that people in modem societies now prefer to forge links or collaborate with those who 
have complementary productive skills and that socioeconomic networks are becoming increas¬ 
ingly heterophilous, but no direct evidence of this has been documented. The availability of big 
data recorded from massively multiplayer online role-playing games (MMORPGs) enables us 
to test social and economic hypotheses and theories—such as this one—in large-scale virtual 
populations (22) and gain a deeper understanding of our social and economic behaviors ( 23 - 31 ). 

We study the evolving collaboration networks of 124 virtual worlds of a popular MMORPG 
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in which the avatars (virtual people) belong to three different professions (warrior, mage and 
priest) with different skills. In a virtual world, each avatar can maintain a list of collaborators 
and the closeness of each pair of collaborators is quantitatively measured by their mutual inti¬ 
macy. More details about the data are found in Material and Methods and SI text. We find that 
avatars prefer to collaborate with others of different professions in these virtual worlds. We then 
propose an economic model to quantify the heterophily by assuming that individuals in socioe¬ 
conomic systems choose collaborators that are of maximum utility. Model calibration shows 
that the evolving curves preference coefficient 7 y(£) that quantifies the degree of collaboration 
preference from /-avatars to 7 -avatars are well separated and illustrate remarkable heterophily: 
7 ij(t) > 1 for avatars belonging to two different professions (/ 7 ^ j ) and 777 (f) < 1 for avatars 
of the same profession (i = j). We further construct two measures to quantify the degree of 
complementarity of virtual societies. We find that social complementarity positively correlates 
with economic output. 


Results 

Empirical analysis 

Fig. Q]4 shows the collaboration networks on day t = 15 of a group of 27 avatars randomly 
chosen from a virtual society filtered by three intimacy thresholds I c = 0, 100, and 2000. There 
are 12 warriors, 5 priests, and 10 mages. If /-avatars are homophilous (neutral, heterophilous) 
in their collaboration-forging process, the proportion of links between /-avatars is greater than 
(equal to, less than) the square of the proportion of /-avatars (0.1975 for warriors, 0.0343 for 
priests, and 0.1372 for mages). For I c = 0, there are 77 links including 15 intra-warrior links, 
4 intra-priest links, and 4 intra-mage links. The proportions of intra-profession links are 0.1948 
for warriors, 0.0519 for priests, and 0.0519 for mages. For I c = 100, there are 48 links in¬ 
cluding 8 intra-warrior links, 1 intra-priest link, and 1 intra-mage link. The proportions of 
intra-profession links are 0.1667 for warriors, 0.0208 for priests, and 0.0208 for mages. For 
I c = 2000, there are 15 links including only one intra-warrior link and no intra-priest and intra¬ 
mage links. The proportions of intra-profession links are 0.0667 for warriors and 0 for priests 
and mages. Hence, the avatars in Fig. Q]4 are heterophious except for priests when I c = 0. 
We will show below that heterophily is not a specific characteristic for these 27 avatars but a 
universal feature present in all the virtual societies. 

Consider a society s on day t whose size is the number N s of s-avatars. The number of 
(s, /)-avatars is denoted N SJ: , where % — 1, 2, and 3 stand for the three professions. Hence 
N s = X/*=i N s.i- The ratio of /-avatars in society s is 

w s ,i = N Sti /N s . (1) 

The average number of j -collaborators of an (s, /)-avatar is f SJJ . Hence, the average number 
of collaborators that an (s,i )-avatar has is f s , t = Y^ =1 fs,ij- The average proportion of j- 
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collaborators in all collaborators of an z-avatar is 



( 2 ) 


Note that q s ^ is the homophily index (15,32). If z-avatars have zero preference for collaborating 
with j-avatars, we have q Sji j = w s j. If z-avatars prefer to collaborate with j-avatars, we have 
q Sji j > w s j. In this case the z-avatars are homophilous when j = i and the -/'-avatars are 
heterophilous when j^z. Fig. QJ? shows that when t — 15 and I c = 100 most virtual societies 
have q s , t j > w s j when i ^ j, but q si j < w s j when i = j. Such heterophilous patterns are 
observed for other values of t and I c as well. 

Similar to the inbreeding homophily index (15, 32), we define the collaboration preference 
index to be 



( 3 ) 


Note that P s/ll is the inbreeding homophily value (15,32). If /'-avatars have no preference to col¬ 
laborate with j-avatars, we have P Sji j = 0. If z-avatars prefer to collaborate with j-avatars, we 
have P S jj > 0. In the latter case, the z-avatars are homophilous when j — i and heterophilous 
when jy^i. Empirical results show that for most virtual societies 1 :> SIJ > 0 when z =4 j, but 
P s i j < 0 when i — j (Fig.[l]Cj. Thus in socioeconomic networks the avatars are heterophilous. 

Fig. HP shows the evolution of preference values P^ averaged over all societies on the same 
day for I c = 100. Although these curves exhibit mild trends, it is evident that the heterophilous 
feature is persistent as the virtual societies develop. 

Quantifying collaboration preference 

Viewing virtual worlds as analogous to real societies, an avatar seeks collaborators when she 
finds it difficult to complete a task alone. The choice of collaborators has a significant influence 
on the completion of the task and it is better to have collaborators with complementary skills. 
Hence the number and skill configuration (or distribution) of an avatar’s collaborators are the 
main determinants of her utility. If the skill configuration in the collaborator list of an avatar is 
optimal, her utility reaches its maximum. If the skill configuration deviates from that optimal 
value, her utility is reduced. Taking into account the fact that maintaining a social network 
incurs a cost, we propose an economic model for the decision-making process of an avatar to 
determine the optimal number of collaborators needed (see Materials and Methods). The model 
essentially maximizes the function 
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where 7 ^ is the preference coefficient. If z-avatars prefer to collaborate with j-avatars, 7 ^ > 1; 
if z-avatars have no preference regarding j-avatars, 7 ^ = 1 ; if z-avatars prefer not to collaborate 
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with 7 -avatars, 7 ^ < 1. For {i,j, k) = {1, 2, 3}, if 7 ^ > 7 ^, then i-avatars prefer 7 -avatars 
over fc-avatars. 

To calibrate the model, we follow and further develop an econometric method presented in 
Ref. (72) (see Materials and Methods). We obtain the values of 7 ^ for each intimacy threshold 
I c on each day t. Fig. [2]4 shows the evolution of preference coefficients 7 ^ for socioeconomic 
networks using the intimacy threshold I c = 100, and Fig. [2^1 shows the average preference 
coefficients over all days. The F-tcsts presented in Materials and Methods show that all the 
results are significant at the 0 . 1 % level. 

All the estimated values of the 7 ^ coefficients are less than 1, while all the 7 ^ values for 
i 7 ^ j are greater than 1. This indicates that the avatars are not seeking same-profession avatars 
but different-profession avatars and are thus heterophilous. In most cases, especially when the 
intimacy threshold I c is not large, the 7 ^ ( I c . t ) values do not have a trend along the evolution 
of virtual worlds. When I c is large, however, we observe an increasing trend in 7 13 (J c ,f) for 
I c = 1000 and 2000, in 723 (/ c , t ) for I c = 1000 and 2000, and in 732 (/ c , t) for h = 500, 1000 
and 2000. We find that the preference coefficients might change with the increase of I c . For 
warriors, 711 and 713 decreases, while 712 increases. For priests, 721 increases, 722 does not 
exhibit evident trend, while 723 decreases. For warriors, 731 increases, 732 decreases, while 712 
increases for large I c values. 

There are also intriguing patterns of relative collaboration preference as quantified by 7 ^ — 
7 ik where i, j and k correspond to the three professions (Fig. [2]S). On average, warriors prefer 
priests over mages and this relative preference enhances when I c becomes greater but reduces 
slightly when t increases for large I c values. Priests prefer mages over warriors when I c values 
are small and prefer warriors over mages when I c values are large. For large J c , priests’ relative 
preference on warriors over mages decreases along time t. Mages prefer priests over warriors 
when I c is small and prefers warriors over priests when I c is large. For large I c , mages’ relative 
preference on warriors over priests also decreases along time t. 


Group complementarity and economic output 

To measure the economic implications of heterophilous preference in socioeconomic networks, 
we investigate the relationship between complementarity of professions and economic perfor¬ 
mance. Consider the socioeconomic network J\f s (I c , t ) of a virtual society with intimacy thresh¬ 
old I c on day t. Economic production utilizes virtual money and goods that are converted to a 
standardized currency (see Materials and Methods). For each member avatar a in J\f s (I c , t), we 
calculate her production output in the week from t — 6 to t, denoted as The economic 

performance of the avatars in JV S (I C , t ) is defined as the output per capita, 


Y s (I c ,t) 
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£ ^ 

a€Af s (Ic,t) 


( 5 ) 
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One measure of profession complementarity can be defined as the sum of preference measures 
between the three types of avatars, 


PsilcJ) = Y.Ps^t). 


( 6 ) 


Alternatively, we can measure complementarity by determining how much the real collaborator 
configuration q s ^ deviates from the optimal collaborator configuration 7 ijW s j (see Materials 
and Methods). The lower the deviation, the higher the degree of complementarity. Thus, we 
have 
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To make these results comparable for different virtual worlds, we investigate the rela¬ 
tive quantities between two societies in the same world, lg(P 2 fc-i/ Pik), lg(C' 2 /, -1 /G' 2 /,) and 
lg(^ 2 fc-i/^ 2 fc). rather than focusing on each society separately. Both measures of complemen¬ 
tarity correlate strongly with the relative economic output when t and / are not large (Fig. [3}4- 
F). For the first few days (small t), most avatars strive to achieve higher levels by implementing 
specific tasks with small economic outputs. Other avatars attempt to obtain high intimacy lev¬ 
els by killing monsters in locations unrelated to economic outputs. In both cases the avatars 
intend to form complementary collaboration networks, but their activities are not focused on 
economic outputs. With the development of a virtual world, the number of active avatars in¬ 
creases and reaches a maximum at time t max and then decays (Fig. EtT). When the activity level 
of a virtual world decreases, the intent of the avatars moves away from production and the col¬ 
laboration structure is increasingly unrelated to economic activities. This is consistent with the 
fact that the spectrum of f max has a distribution similar to the significant correlations between 
complementarity and economic output (Fig. [3]f/). 

Discussion 

Overwhelming empirical evidence has shown that most social networks are homophilous. The 
probability that two nodes will connect is higher if they share similar traits. Our analysis of 
virtual worlds in which division of labor is operative demonstrates the important role of com¬ 
plementarity. In those socioeconomic networks individuals have the motivation to cooperate, 
and in the formation of the network individuals exhibit a heterophilous preference for those 
with complementary productive skills. Although mapping human behavior in virtual worlds 
to real-world human behavior is a subtle process (33), we believe that they share an intrinsic 
commonality because avatars in virtual worlds are, in fact, controlled by real-world people. 
In particular, avatars consciously form teams to accomplish tasks more successfully and ef¬ 
fectively. More generally, growing evidence shows significant similarities in the behaviors of 
online avatars and real-world humans (22, 34-41) 
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In reality, human’s preference is multidimensional in their traits (13). The situation in virtual 
societies is a little different. Indeed, the way people interact with each other has significantly 
changed from the old days, particularly due to the impact of internet. In modem time, people 
can meet through internet in the virtual world instead of physically getting together to dine, 
drink, and talk to forge ties. Personal traits become less important in virtual societies while 
avatars’ profession skill is identified as a dominating trait in virtual society. 

The economic model proposed in this work is different from the one in (12). The essential 
difference is in the assumption of the utility function. Also, we have used a modified method 
of model calibration. Moreover, our model allows us to determine not only if z-avatars are 
homophilous or heterophilous but also the preference of one type of avatars to any other type of 
avatars. Hence, our model is more general and can be applied to other systems. 

The relationship between social networks and economic output has been studied previously. 
It has been found, for example, that the diversity of individual relationships within a community 
strongly correlates with the economic development of the community (42) and is directly asso¬ 
ciated with higher productivity for both individuals and the community (43, 44). Because, to 
date, detailed real data at the population level of societies have been unavailable, this correspon¬ 
dence between professional skill and economic performance has not been quantified. Here we 
have begun to fill this data gap and also to highlight the usefulness of virtual worlds in carrying 
out research in economics and sociology (22). One direct indication of our findings is that if 
a team leader or a firm manager recruits new members according the complementarity of their 
skills, the team’s productivity will increase and the firm’s economic well-being grow. 
Supplementary Material accompanies this paper at http: / /www. scienceadvances . org/. 


Materials and Methods 

Data description 

We use a huge database recorded from K = 124 servers of a popular MMORPG in China to 
uncover the patterns characterizing virtual socioeconomic networks. In a virtual world residing 
in a server there are two opposing camps or societies. Two avatars can choose collaborators, 
and a measure of closeness called intimacy is assigned to the collaboration link. When two 
collaborators in the same society collaborate to accomplish a task, their intimacy level increases. 
Two avatars from different societies can also collaborate, but their intimacy level remains zero. 
Hence the social networks of the two camps are essentially separate. We can regard the two 
camps as two societies, thus giving us S — 248 virtual societies. For convenience, s = 2k — 1 
and s = 2k stand respectively for the two societies in the same virtual world k. Two avatars are 
defined as collaborators if they both are on the collaborator list and their intimacy exceeds I c . 
We consider many temporal collaboration networks. On day t in a virtual society s, a network 
Af s (l ci t ) is a network in which the intimacies of all edges are no less than a threshold I c , which 
can be disconnected (Fig. Q]4). 
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In each society there are three professions (warrior, priest, and mage). We use subscripts 1, 
2, and 3 to stand respectively for the three professions: warrior, priest and mage. For simplicity, 
we define several notations as follows. An s-avatar is an avatar belonging to society s. An 
/-avatar is an avatar having profession /. Similarly, an /-collaborator is a collaborator having 
profession i. An (s, /)-avatar or (.s, /)-collaborator is an /-avatar or /-collaborator in society s. 


An economic model of collaboration formation 


If there is no collaboration preference, the proportion of (s, j)-collaborators that an (s, /)-avatar 
has is identical to the proportion of j-avatars in the group, that is q Sii j = w s j. Hence the 
number of (s, j)-collaborators of an (s, z)-avatar is j ™ 1 = f s ^w s j. In a virtual society with a 
preset division of labor, an avatar prefers collaborators with complementary skills. We assume 
that, for an (s, /)-avatar, there is an optimal configuration of collaborators with different skills, 
f°Tj = where the preference coefficients 7 ^ are independent of society s. The utility 

of an (s, /)-avatar increases when her/his real number f S:lJ of (s, j)-collaborators approaches the 
optimal value and reaches its maximum (7™ ax = bf ^ when her/his collaborator configura¬ 
tion is optimal such that f S:i j = /°?j. According to the law of diminishing marginal utility, we 
have f3 < 1 . Therefore, the utility function of an (7 /)-avatar is 


U s ,i = b s fg,i ~ a a f%XZ. 


(3—a 


( 8 ) 


where 




'y X VsM 7 ij w s,j 

.3 =1 


2(/3-c) 


( 9 ) 


and 7 ij is the preference of (s, /)-avatars for (.7 7) -avatars. If /-avatars do not have any prefer¬ 
ence on j-avatars such that q syj = w s j for all societies, we have 7 y = 1. If /-avatars prefer 
j-avatars, we have 7 ^ > 1. If /-avatars prefer not to collaborate with j-avatars, we have 7 ^ < 1. 
For {/, j, k) = (1, 2, 3}, if 7 y > 7,;/,, then /-avatars prefer j-avatars over /c-avatars. To maintain 
a collaboration network of size f s/l . the (s, /)-avatar suffers a cost proportional to f SJ (12), 


Ms,i = Csf V.. 


( 10 ) 


According to the above model, the overall utility in the decision-making process is 


D 8 ,i(fsJ = U 8ii - M 8ii = ( b s - c 8 )f J 




(3—a 


(ii) 


Model calibration 


An (s,/)-avatar solves the following decision-making problem of how many collaborators to 
have 


ma xD Sji (f s>i ). (12) 
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It follows that 


fa,i = K«/ (b a ~ c a )0] ?-<* X S)i . (13) 

Note that the 7 ^ values are affected only by the professions and remain the same for different 
societies. This enables us to estimate the parameters. 

The solution (fT3l) denotes the average behavior (decision) of all avatars having the same 
profession in a given society. If we consider an arbitrary avatar a , we must add a noise term (12), 

fs,i-,a = [ a s a/(b s - c s )/3 ]^ X Sti + e a , (14) 


which means that the “realized” number of collaborators avatar a has is the sum of a universal 
(or systemic) term and an idiosyncratic error term. The error term is assumed to have mean 0 and 
variance a 2 . Note that this assumption states that the variance of any avatar of any profession is 
the same. 

We denote N s as the size of society s and u;. s/ , as the fraction of /’-avatars in society s. Hence 
the number of /-avatars in society s is N s w s j, and the expectation of the aggregated number of 
collaborators that /-avatars have in society s is N s w S:i f s According to Eq. (fl4l) . we have 


Xs^s,ifs,i X S W S 


a*a 


_(b s - c,)j 8 


Xs,i T E s i, 


where E Sji has mean 0 and variance ( N s w S)i a) 2 = ( N s w sti ) 2 a 2 . 
It follows that, for i ^ j, 

m s,iX s fs,i E s j y Ws,jX s f s j E 


s,3 


^s,iX s X s i 


X S X ggj 


(15) 


(16) 


For society s, the error is 


^s,ij ^s,iX s f s,i^s,jX s X s j W s jX s f s jVJ s jX s X s j 

XgU!gjWgj(fgjX s j fgjXgi) 

Xs. is. j X s X s; j EgjWgjNgXg^. 


(17) 

(18) 


According to Eq. (fT8l) . we find that the mean of is 0 and the variance is 4> Sj ij<J 2 , where 

= (w s ,{N S ) 2 (w 3 jN g X 3:j ) 2 + ( w 3 jN 3 ) 2 (w s:i N s X S:i ) 2 

= /v>, 2 X7(Wi) 2 + W’l (19) 

Thus the normalized variable ^ 2 , tJ /<p S: ij has mean 0 and variance a 2 for any society s. The sum 
of squared errors (Qfj) over all societies in the sample is 


Ql-f 


S =1 



4*s,ij 


{fs,jX s ,j — f s ,jX S) i ) 2 

h {x s ,f+(x s ,f 


( 20 ) 
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which is independent of N s as expected. However, Q 2 j is dependent on w s> i, which is consistent 
with the setup of our model but different from the model in Ref. (72). Thus the total sum of the 
squared errors is 

Q 2 — Q 12 + Q13 + Qh- ( 21 ) 

One can see that a s , b s and c s could be society-specific and are not included in the final objective 
function of model calibration. 

For each pair of I c and t, a society is excluded in model calibration if the number of avatars 
having at least one collaborator is less than 500 to ensure that e a has enough realizations. Chang¬ 
ing this threshold around 500 results in same results. In addition, if the number of societies 
included in a model is less than 50, we do not calibrate the model because the model has 10 
parameters. 

To find the solution to the minimization of Q 2 , the taboo search algorithm is adopted (45). 
The solution space is restricted to 0 < 7 ^ < 2 for i , j e {1, 2, 3} and 2 {p- a ) > 0. Because there 
are 10 free parameters, it is not easy to reach the global minimum. We thus perform a taboo 
search in each cell of a 9-dimensional lattice of size 2 9 with the constraint that 0 < 7 ^ < 1 or 
1 < 7 ^ < 2 . The parameters in certain cell corresponding to the minimum of Q 2 in all cells are 
obtained as the solution. The normality assumption of fitting errors has been verified by QQ- 
plots, which rationalizes the setup of the model. We note that the partitioning of the solution 
space into a 9-dimensional lattice of size 2 9 is very important. If we perform the taboo search 
directly, the resulting Q 2 value is significantly larger and the three preference curves 77 (t) for 
each i are not well separated around 7 ^ = 1 (cf. Fig. [2). 


Significance tests 

To test whether the preference coefficient 7 \j of 7-avatars to j-avatars is significantly different 
from the no-preference case, we perform F-tcsts using the null hypothesis 


H 0 - 7 ij = h (j'e {1,2,3}. 


( 22 ) 


Following Ref. (72), the F-statistic is 


F 


SSR 

P 


con S S 7? U ncon 


uncon Peon 


SSR uncon 


Puncon 


(23) 


where SSR is the sum of squared residuals of the best-fit calibration, p is the number of model 
parameters, n is the number of observations, while the subscript “con” indicates the constrained 
model under the null hypothesis and the subscript “uncon” the unconstrained model. 


Economic output of individuals 

There are two virtual currencies, Xingbi and Jinbi. Xingbi cannot be produced by an avatar’s 
activity and can only be bought from the system, which has an approximately stable exchange 
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rate in reference to the Chinese currency Renminbi. Xingbi is thus a universal currency across 
different virtual worlds. Jinbi, on the other hand, is produced by the economic activities of the 
avatars. There is a built-in exchange platform in each virtual world so that avatars can exchange 
Xingbi and Jinbi. In this way, there is a real-time exchange rate from Jinbi to Xingbi. 

An avatar can produce virtual items (e.g., weapons, clothes, and medicines) and a limited 
amount of the virtual currency Jinbi. We convert the produced items and Jinbi to Xingbi to 
obtain the real economic output of each avatar on each day. There is a marketplace in each 
virtual world in which avatars can sell their items that are priced in Xingbi or Jinbi. The price 
of an item is determined by the average price of all the trades in the marketplace on a given day. 
Each produced item can thus be measured in Xingbi. 
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Figure 1: Empirical evidence of heterophily in the socioeconomic networks of virtual societies 
on a typical day t = 15. Warriors, priests and mages are marked respectively in cyan, red 
and blue. (A) Networks of 27 avatars randomly chosen from a virtual society filtered by three 
intimacy thresholds I c = 0, 100 and 2000 (top to bottom). ( B ) Dependence of q s , tJ on relative 
size w s j for all virtual societies for I c = 100. In each plot, there are three well isolated clusters. 
For most societies, q s , t j > w s j when i ^ j and q, %lJ < w S: j when i = j. (C) Dependence of 
preference measure P S)i j on relative size w s j for all societies for I c = 100. There are also three 
well separated clusters in each plot. For most societies, P S}i j > 0 when i ^ j and P S/IJ < 0 
when i = j. ( D) Evolution of the averaged preference measure P K /, 3 over all virtual societies 
for I c = 100. The preference measures are roughly persistent. 
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Figure 2: Preference coefficients 7 ^ for socioeconomic networks with the intimacy threshold 
being I c = 100. (A) Daily evolution of the nine preference coefficients 7 ^ with i,j e {1, 2, 3}. 
The color of a point (t, 7 ^) is determined by j: cyan, red and blue for j = 1,2 and 3, respec¬ 
tively. The nine points for a given t were determined simultaneously in one calibration. ( B ) Box 
plots of 7 ij shown in A. 
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Figure 3: Relation between complementarity of collaboration network and economic output. 
(A) Examples of correlations between lg(P 2fc _i/P 2fc ) and lg(Yy_i /Y-ik)- ( B ) Examples of cor¬ 
relations between \g(C 2 k-i/C 2 k) and lg(Y 2fc _i/Yy). (C) The p-value of the correlation between 

1 g (/2 a- _i /P 2 k ) and lg(r 2fc _i/y 2 ,) for different values of I c and t. A give grid (t, I c ) is colored 

as red or yellow if the correlation is significant at the 0.001 level or the 0.01 level. Otherwise, 
the grid is colored as green. ( D) The />-value of the correlation between lg(C 2fc _i/C 2fc ) and 
lg(^ 2 fc—i/^ 2 fc)- (E) Correlation coefficient p between lg(P 2 fc_i/P 2 fc) and lg(Y 2k -i/Y 2k ) for dif¬ 
ferent values of I c and t. The correlation coefficient is set to be zero is the correlation is insignif¬ 
icant at the 0.01 level. ( F ) Correlation coefficient p between lg(C' 2 /,._i /Gy) and lg(y 2fc _i/y 2fc ). 
(G) Evolution of the number of active avatars in different virtual worlds. ( H ) Histogram of t max 
which is the date that a virtual world has historically the maximum active avatars. 


17 



























































